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Abstract 

When lepton flavour effects in thermal leptogenesis are active, they introduce important dif- 
ferences with respect to the case in which they are neglected, the so-called one-flavour ap- 
proximation. We investigate analytically and numerically the transition from the one-flavour 
to the two-flavour case when the r-lepton flavour becomes distinguishable from the other two 
flavours. We study the impact of the oscillations of the asymmetries in lepton flavour space 
on the final lepton asymmetries, for the hierarchical right-handed neutrino mass spectrum. 
Flavour oscillations project the lepton state on the flavour basis very efficiently. We conclude 
that flavour effects are relevant typically for Mi < 10^^ GeV, where Mi is the mass of the 
lightest right-handed neutrino. 
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1 Introduction 



Baryogenesis through Leptogenesis [T] is a simple mechanism to explain the observed baryon 
asymmetry of the Universe. A lepton asymmetry is dynamically generated and then converted 
into a baryon asymmetry due to {B + L)-violating sphaleron interactions [2] which exist within 
the Standard Model (SM). A simple scheme in which this mechanism can be implemented 
is the 'see-saw' (type I) model of neutrino mass generation jj]- In its minimal version it 
includes the SM plus two or three right-handed (RH) heavy Majorana neutrinos. Thermal 
leptogenesis [H El E] can take place, for instance, in the case of hierarchical spectrum of 
the heavy RH Majorana neutrinos. The lightest of the RH Majorana neutrinos is produced 
by thermal scattering after inflation. It subsequently decays out-of-equilibrium in a lepton 
number and Charge and Parity (CP) violating way, thus satisfying Sakharov's conditions [7j. 
On the other hand, the see-saw mechanism of neutrino mass generation [3j provides a natural 
explanation of the smallness of neutrino masses: integrating out the heavy RH Majorana 
neutrinos generates a mass term of Majorana type for the left-handed flavour neutrinos, which 
is inversely proportional to the large mass of the RH ones. 

The importance of the lepton flavour effects in thermal leptogenesis has been recently 
realized in [HI [9l [ini E] . The dynamics of leptogenesis was usually addressed within the 'one- 
flavour' approximation. In the latter, the Boltzmann equations are written for the abundance 
of the lightest RH Majorana neutrino, A^i, responsible for the out of equilibrium and CP- 
asymmetric decays, and for the total lepton charge asymmetry. However, this 'one-flavour' 
approximation is rigorously correct only when the interactions mediated by charged lepton 
Yukawa couplings are out of equilibrium. Assuming for the moment that leptogenesis takes 
place at temperatures T ~ Mi, where Mi is the mass of A'^i, and that the RH spectrum is hi- 
erarchical, the 'one-flavour' approximation holds only for T ~ Mi ^ 10^^ GeV. For Mi > 10^^ 
GeV, all lepton flavours are not distinguishable. The lepton asymmetry generated in A'^i decays 
is effectively 'stored' in one lepton flavour. However, for T ~ Mi ^ 10^^ GeV, the interactions 
mediated by the r-lepton Yukawa couplings come into equilibrium, followed by those medi- 
ated by the muon Yukawa couplings at T ~ Mi < 10^ GeV, and the notion of lepton flavour 
becomes physical. Flavour effects are important because leptogenesis is a dynamical process, 
involving the production and destruction of the heavy RH Majorana neutrinos, and of a lepton 
asymmetry that is distributed among distinguishable flavours. Contrary to what is generically 
assumed in the one-flavour approximation, the AL = 1 inverse decay processes which wash 
out the net lepton number are flavour dependent. The asymmetries in each lepton flavour, are 
therefore washed out differently, and will appear with different weights in the final formula 
for the baryon asymmetry. This is physically inequivalent to the treatment of washout in the 



1 



one-flavour approximation, where the flavours are taken indistinguishable. 

The impact of flavour in thermal leptogenesis has been recently investigated in detail in 
O Uni [m [12l [13] , including the quantum oscillations/correlations of the asymmetries in lepton 
flavour space [U]. The interactions related to the charged Yukawa couplings enter in the 
dynamics by inducing nonvanishing quantum oscillations among the lepton asymmetries in 
flavour space [9j. Therefore the lepton asymmetries must be represented as a matrix Y in 
flavour space, the diagonal elements are the flavour asymmetries, and the off-diagonals encode 
the quantum correlations. The off-diagonals should decay away when the charged Yukawa 
couplings mediate very fast processes. The Boltzmann equations therefore contain new terms 
encoding all the information about the action of the decoherent plasma onto the coherence 
of the flavour oscillations: if the damping rate is large, the quantum correlations among the 
flavours asymmetries are quickly damped away. If leptogenesis takes place when the charged 
Yukawa couplings do not mediate processes in thermal equilibrium, the quantum correlators 
play a crucial role to recover the one- flavour approximation. On the other hand, if leptogenesis 
takes place when the charged Yukawa couplings mediate processes well in thermal equilibrium 
quantum correlations play no role in the dynamics of leptogenesis. 

The goal of this paper is to study the transition from the one-flavour to the two-flavour case. 
In the case of hierarchical RH mass spectrum, the baryon asymmetry is directly proportional 
to the mass Mi of the lightest RH neutrino. A large enough baryon asymmetry is obtained 
only for a sufficiently large value of Mi. Therefore, we will restrict ourselves to the transition 
from the one-flavour state, to be identified with the total lepton number, to the two-flavour 
states, to be identified with the r lepton doublet £r and a linear combination of the fj, and 
e doublets. The most interesting region is for values of masses of the lightest RH neutrino 
centered around Mi ~ 10^^ GeV where we expect the quantum correlators to play a significant 
role in projecting the lepton state on the flavour basis and, eventually, in the generation of the 
baryon asymmetry. Studying the details of the transition is relevant to understand if it is a 
good approximation to compute the baryon asymmetry just solving the Boltzmann equations 
with only the diagonal entries of the matrix Y for the lepton asymmetries (as usually done 
in the recent literature for the flavoured leptogenesis [IH [T2t I13j ) and neglecting altogether 
the off-diagonal entries. We would like to see under which conditions on the leptogenesis 
parameters the full two-flavour regime is attained. 

The paper is organized as follows. In Section [2] we summarize the general framework and 
the Boltzmann equations. In Section [3] we describe in detail the one-flavour limit, while the 
two-flavour limit is described in Section Section [5] contains the main body of our results; we 
present both analytical and numerical results for the various regimes. Finally our conclusions 
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are contained in Section [6] together with some comments. 

2 Two-flavour Boltzmann equations 

The lagrangian we consider consists of the SM one plus three RH neutrinos Ni (i = 1,2,3), 
with Majorana masses Mj. Such RH neutrinos are assumed to be heavy (i.e. with masses weh 
above the weak scale) and hierarchical (Mi <C ^2^3), so that we can safely focus our attention 
on the dynamics of A'^i only. The interactions among RH neutrinos, Higgs doublets H, lepton 
doublets ia and singlets Cq, (a = e, /i, r) are described by the lagrangian 

^int = KaNii^H + Ke^i^H' + ]^M,N,Ni + h.c. , (1) 

with summation over repeated indeces. The lagrangian is written in the mass eigenstate 
basis of RH neutrinos and charged leptons. The interactions mediated by the charged lepton 
Yukawa couplings are out of equilibrium for T ~ Mi ^ lO^'^ GeV. In this regime, flavours are 
indistinguishable and one can perform a rotation in flavour space to store all the asymmetry in 
a single flavour. At smaller temperatures, though, this operation is not possible. The r flavour 
becomes distinguishable for T ~ Mi < 10^^ GeV. As we already discussed in the Introduction, 
we will restrict ourselves to the study of the transition occuring around T ~ Mi ~ 10^^ 
GeV. This choice is motived by the following considerations. In the case of hierarchical RH 
mass spectrum, the baryon asymmetry is directly proportional to the mass Mi of the lightest 
RH neutrino. Therefore, a large enough baryon asymmetry is obtained only for a sufficiently 
large value of Mi. Since the transition which makes the /i flavour distinguishable occurs at 
T ~ Ml ~ 10^ GeV, the corresponding value of Mi is generically too small to provide a baryon 
asymmetry in the observed range. Therefore, we will study the transition from the one-flavour 
state, to be identified with the total lepton number stored in the lepton doublets, to the two- 
flavour states, to be identified with the r lepton doublet ir and a linear combination of the 
and e doublets (which at temperatures between 10^ and 10^^ GeV are indistinguishable), 

1 /2 

h = (Ale4 + \lM I (lAlel' + \\l/) . 

Having therefore in mind the transition between a one-flavour and a two-flavour system, we 
study a toy model with two lepton doublets a = 1, 2 and generically represent the lepton asym- 
metry matrix by a 2 x 2 density matrix Y given by the difference of the density matrices for the 
lepton and anti-lepton number densities (normalized to the entropy density s). The diagonal 
elements are the lepton asymmetries stored in each flavour while the off-diagonal elements 
describe the quantum correlations between different flavours. The total lepton asymmetry is 
given by the trace of this matrix. 
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{id + 7AL=l) ( - 1 ) e«/3 - l^Si{lD + 7AL=1, i^} 



a/3 



In order to follow the evolution of the lepton asymmetry, one needs to write down the 
equations of motion for the matrix Y . The proper evolution equations for the matrix Y has 
been found and discussed in [9], neglecting the transformations to bring the asymmetries in 
the lepton doublets to the the SM conserved charges = {B/3 — L^), where is the total 
lepton number in a single flavour. Including these transformations only change the final result 
by a factor of order unity and therefore we will also neglect them for the sake of presentation. 
The interactions mediated by the Yukawa couplings are also taken into account. We will 
assume a large hierarchy between the Yukawa couplings (which holds for the realistic case, 
since hr ^» ^^i.e)- 

The system of Boltzmann equations for the generic components Y^/s of the density matrix, 
as a function of the variable z = Mi/T, read[f| 

dY^^ 1 

- cT2Re(A) + ai|Im(A)|] Y^^ , Y^^ = Y;^ , (2) 
while the Boltzmann equation for the A^i abundance (Yn^) is 

where the equilibrium A'^i abundance is given by Y^{z) = ^2^X2(2:), and g^, is the number of 
effective degrees of freedom in the thermal bath. Notice that we have included the contribution 
to the CP asymmetry from the AL = 1 scatterings [11] . 

We remark that to obtain Eq. ([2|) we have assumed that the lepton asymmetries oscillate 
with an approximately momentum-independent frequency. The oscillation frequency in flavour 
space depends on the energy (momentum) of the leptons and, within one oscillation timescale, 
leptons are involved in many momemtum-changing interactions caused by the fast, but flavour- 
blind, gauge interactions. Our assumption amounts to adopting the thermally averaged energy 
{E) to estimate the oscillation frequency. In other words, we have approximated the integral 
J iEdt with i{E) J dt along the path from one lepton number violating interaction to the next. 
This approximation is well justified in [T3], where it has been shown that fast gauge interactions 
do not affect the coherence of the flavour oscillations. 

Before discussing the Eqs. ([2]) and (l3|), we explain the various quantities appearing in them. 
The matrix (70)0/3 represents the thermally averaged A^i-decay rates ant it is given by 

(7z))a/3 = 7D7YTt^ = 7D^^-n— [2 > (4) 



^As usual, {, } stands for anti-commutator while the cr's are Pauli matrices. 
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normalized in such a way that the total decay rate is the trace of the matrix. The AL = 1 
scatterings were also included in the equations (see [11] for a discussion about this point). The 
thermally averaged interaction rate matrix (7al=i)«/3 has the same form as (713)0-/3 in (j3|) with 
7£) replaced by the total scattering rate 7al=i- The explicit expressions for the total rates 
and 7AL=i can be found in the literature (see e.g. [1]). 

It is possible to generalize the usual decay parameter to the two-flavour case. The natural 
definition is a 2 x 2 matrix 



K, 



Ta/J 



(5) 



z=l 

)af3 
^ycq Ki(2) ' 



H 

where 

and Kj(z) are modified Bessel function of the second kind. The trace of K^^p will be denoted 
by K = J2a^c,a- 

The CP-asymmetry matrix is given by [S]: 



'"'^ 16vr 

where the loop function / is [T5 



(7) 



m 

Notice that 



l-(l + x)log( 1 + - ) + ^ 



1 - X 



(8) 



^af3 — ^I3a (9) 



and the normalization is such that the trace of the CP asymmetries reproduces the total 
CP asymmetry produced by the decays of the lightest RH neutrino A^i, in the single-flavour 
approximation 



(10) 

\ " / \/vi-i 

If m denotes the heaviest light neutrino mass (= matm for the non-degenerate case) then the 
entries of the CP-asymmetry matrix are subject to the bounds [9] 



where v is the vacuum expectation value of the Higgs doublet. 
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The A parameter accounts for interactions mediated by the dominant Yukawa couphng, 
which from now on we denote by hi. It is given by 



^ uji- iTi 



(12) 

T=Mi 



H{Mi) 

having defined the thermal mass oji ~ h\T /\& and the interaction rate Fi ~ 8 x 10~^/ifT |16| . 
The dependence on Mi is easily made explicit: 

Re(A) ~ 4 X , Im(A) ~ -5 x W'^hj^ , Re(A) ~ 10 |Im(A)| , (13) 

Ml Ml 

where Mp = 1.2 x 10^*^ GeV is the Planck mass. In the realistic case, one should identify hi 
with hr- The flavour 1 will therefore become distinguishable when Mi ^ 10^^(/ii//it-)^ GeV. 

The parameter A will play a crucial role in what follows. It contains all the informa- 
tions about the action of the decoherent plasma onto the coherence of the flavour oscillations. 
Changing the parameter A, that is changing the value of the mass Mi, and assuming that 
leptogenesis takes place at a temperature T ~ Mi, one can analyze the various regimes: for 
|A| <C 1, the Yukawa coupling hi does not mediate processes in thermal equilibrium and one 
expects therefore that the one-flavour approximation holds. In this regime the off-diagonal en- 
tries Ya/3 are expected to be nonvanishing. For |A| ~ 1 the transition between the one-flavour 
and the two-flavour states takes place. For |A| ^ 1 the transition is occured, there are two 
flavours in the system and one expects the off-diagonal entries in the matrix Y to be decaying 
very fast since the quantum correlations among the flavours is efficiently damped away by the 
decoherent interactions with the plasma. 

It is simpler to work with the Boltzmann equations obtained from ©-([3I) by eliminating 
the thermally averaged rates in favor of the decay parameter matrix K^j^p and two functions, 
fi{z) and f2iz), which account for the AL = 1 scatterings in the Ni thermalization and in the 
wash-out of the asymmetry, respectively (see [Hill]). Their asymptotic behaviours are 

/.(^)-{ (14) 

for z < 1 , 




and 



. / N I 1 for z > 1 

f2iz)^< (15) 



where g^z^ ~ 0.1 parametrizes the strength of the AL = 1 scatterings and uk = 4/3 (2) for 
the weak (strong) wash out case. A good approximation to the total wash-out term (inverse 
decays and AL = 1 scatterings) at small z is given by ~ l{)^^aKK . 
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After a short manipulation the Boltzmann equations read 

1 



YLp = -y^,e„^--/i(z){K,y},^- [a2Re(A) + ai|Im(A)|Jy,^, (16) 



= -zK-p^h{z){YN,-Y;^^), (17) 



- Ki(^) 
K2(^)- 

where primes denote derivatives with respect to z and h{z) = ^z^Ki(z)/2(z). These equations 
are the starting point of our analysis. Although they are just classical equations, they reproduce 
the correct expected limits (as shown in the next two Sections) and also have the virtue of 
providing information on the transition between the one-flavour and the two-flavour regimes. 

3 The one-flavour limit 

In this section we deal with the one- flavour limit, corresponding to |A| ^ 1. More precisely, 
inspecting Eq. (jl6p . one learns that the quantum correlators need to be accounted for ifl 

\K\<.\h{z)K. (18) 

which implies 

Ml 



1012 GeV; ^ i^/i(z) ■ ^^^^ 
This condition has to be satisfied at the time when the asymmetry is generated. In the 
weak wash-out regime, K "^1, and supposing that the initial abundance of RH neutrinos is 
vanishing, the production of the baryon asymmetry takes place at some z^l. Since the 
wash-out term for K "^1 is always smaller than unity, we conclude that in the weak wash-out 
regime the one-flavour limit is reached for Mi > IO12 GeV. 

In the strong wash-out regime, A' ^ 1, the baryon asymmetry is generated at some z ~ 
Ixi K + (5/2) In z > 1 when Kh{z)/2 ~ 1. Since the wash-out function Kh{z)/2 is larger than 
unity for z z, we conclude that in the strong-wash out regime the condition (jlSp implies 
|A| < {l/2)Kh{z) ~ 1, that is Mi > IO12 GeV. 

Under the conditions that the A-terms may be dropped in Eq. ()16p . the latter reads 



Yaa = -y'Ni^aa " [Ka/BYpa + KpcXap] " h{z)KaaYaa , (20) 

YLp = -Y^,e^p-]^h{z)Tv(X)K^p-hi{z)KY^f,, (21) 



^ We thank P. Di Bari for sharing with us prior to publication his paper in collaboration with Blanchet and Raffelt 
|17j where similar considerations have been presented. 
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with a ^ 13 and no summation over repeated indices. Notice that these equations are imphcit, 
since the trace of Y appears in the right hand side. Now, we perform an ad hoc rotation in the 
flavour space. The quantities referred to the new basis will be denoted by a 'hat'. In general, 
we are free to rotate the lepton doublets by a unitary matrix A: 

L = (22) 

{AA^ = 1) and this is equivalent to a basis change in the flavour space. A useful choice for A 
is 

A = — ^ ( ] , (23) 

\/[AAV V -(^12)* (All)* J 

where [AA^Jn = |Aiip + |Aiip = [AA^Jn by the unitarity of A, which leads to the rotated 
Yukawa couplings: 

X = + ° V (24) 

^A^^^^l \ (MY hi + (Ai2)'A22 det[A] J ' 

with det[A] = A11A22 — Ai2A22- The zero entry makes manifest that A''i is coupled only to 

The matrices Kaf^ and e^/? in the new basis are obtained by replacing A — > A; in particular, 
one finds 

Ku=K, = K21 = K22 = (25) 

en = ei , £22 = . (26) 



Thanks to these relations, the equations for the diagonal components (j20j) give 122 = 0, so 

the lepton asymmetry is concentrated on the lepton ii only and it evolves according to the 
equation 

l^^ = -y^ei-/i(z)ify^^, (27) 

which exactly reproduces the Boltzmann equation for the one single flavour. The latter can be 
identified with the total lepton asymmetry, that is with the trace of the lepton asymmetries. 
The total lepton asymmetry in the lepton doublets is indeed the only quantity which treats 
indistinguishably all the flavours. 

4 The two-flavour limit 

Let us now turn to the opposite regime where the A terms are important, i.e. we are in the 
full two-flavour regime. Again, we split (I16p in equations for the diagonal and off-diagonal 
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components of Y: 

YL = -Yl^,eaa-^h{z)[K^pY(^a + K(^aYap]-h{z)K^aYa^, (28) 



Ya/s = -Y^,e^p--h{z)TT{Y)K^p 



h{z)K+{\l^{k)\ + {<y2) 



Yo^p , (29) 



with a ^ (3 and no summation over repeated indices. The A terms appear in the wash-out 
of the off-diagonal elements. Therefore, the solutions of (|29p will contain exponential factors 
of the form e*^^. The real part of A leads to oscillating behaviours, while the imaginary part 
controls the damping. The latter is originated by the decoherence effect of the high temperature 
plasma on the flavour oscillations: if Yukawa coupling hi mediates processes which are fast 
enough, the correlations between different flavours are rapidly lost. Such correlations are 
encoded in the off-diagonal components of the lepton asymmetry density matrix Y . As long 
as the off-diagonal entries become negligibly small, Eq. (i28|) reduces to that studied in [llj . 
where the flavours are considered as completely decoupled and the system of equations reduces 
to two equations for the diagonal entries of the Y matrix. More in detail, we can say that the 
two-flavour state is reached when the oscillations are efficiently damped, i.e when the following 
condition holds 

|Im(A)| > \h{z)K (30) 

or 

' Ml \ ^ 2 



1012 GeV ; ~ Kh{z) ' ^^^^ 
around the point when the baryon asymmetry in a given flavour a is generated. In the weak 
wash-out regime for all flavours, Kaa,K < 1, the flavour asymmetry is generated at ^ 1 
and the function {l/2)Kh{z) is always smaller than unity. Therefore, we obtain that the two 
flavour regime is dynamically relevant for Mi ^ 10^2 GeV. 

In the strong wash-out regime for all flavours, K^a^K ^> 1, the condition (j3ip on the 
mass of the RH neutrino is Mi < (Kaa/K) 10^2 GeV for ~ In K^a + (5/2) In z^ > 1. The 
most stringent bound is obtained for the smallest K^a, which corresponds to the smallest 
wash-out. Of course the bound should be applied only if the same flavour gives also the 
largest asymmetry. This depends upon the CP asymmetry eaa- In particular, if K^a takes 
the smallest value compatible with the strong wash-out, ~ 3 and if the CP asymmetry 
eaa is the largest, then one obtains the most stringent bound. Mi < {3/K)10^'^ GeV. 

In the case of strong wash-out for some flavour a, Kaa ^ 1, but weak wash-out for some 
other flavour /3, Kjs/s ^ 1, the asymmetry in the flavour /3 is generated at z = 0{5) [TT] and the 
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condition on the mass of the hghtest RH neutrino is given by Mi < {10 / K)10^'^ GeV, provided 
that the final baryon asymmetry is mainly generated by the flavour /3. If this is not the case, 
one should apply the condition Mi < {K^a/K) 10^^ QeV ~ 10^^ QgV. 

Let us close this section with a comment. We expect the bounds obtained in this section 
comparing rates to be in fact too restrictive. They have been derived just comparing the 
rate of the AL = 1 inverse decays and scatterings with the rate of damping of the flavour 
oscillations. However, the real dynamics is more involved. For instance, the flavour oscillations 
are characterized by a rapidly oscillating behaviour. The oscillation rate is dictated by |Re(A)| 
which is a factor about ten larger than the damping rate of the flavour oscillations. Re (A) ~ 
10| Im(A)|. This is relevant because computing the flavour asymmetries involves integrals over 
time. Since the flavour oscillations decay and also have an oscillatory behaviour, this restricts 
the range of time integration, thus leading to a suppression of the contribution from the flavour 
oscillations. We therefore expect the influence of the the flavour oscillations to disappear even 
in the vicinity of Mi ~ 10^^ GeV. Our numerical results support this expectation. 

5 The transition between the one- and the two-flavour 

case 

Having elaborated about the two extreme regimes, we now investigate what happens in the 
intermediate region where the one flavour - two flavours transition takes place. To achieve this, 
we perform an analytical study of the solutions of (j28p and (j29p , in two representative regimes 
of -RT's, showing also some numerical simulations to enforce our findings. In the figures we 
will present two different quantities which may serve as indicators of the transition. The first 
quantity is Ya^/ ^aa)dcc which is the ratio between the flavour asymmetry Y^a in the flavour 
a computed solving the full system of Boltzmann equations (j28p and ()29p over the same 
asymmetry {Yaa)Aec computed neglecting the off-diagonal terms in the same equations. This 
ratio should tend to unity in the full two-flavour regime because the off diagonal correlators 
have been efficiently damped out. The second indicator is the ratio of the the trace of the 
2x2 matrix Y , Tr[y] computed solving Eqs. (j28p and (|29p and the asymmetry computed 
in the one-flavour approximation, li-flavoun assuming a single flavour with CP asymmetry ei 
and wash-out parameter K = Kn + K22- This ratio should tend to unity in the one- flavour 
regime, when the off-diagonal terms are not damped. 
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5.1 Strong wash-out regime for all flavours 



In this case Kii,K22 ^ 1- This imphes that the A'^i abundance closely follow the equilibrium 
abundance, ~ (^jVi)' ~ ~'^^{^)/^- 'T^^ integrals giving the lepton asymmetries are 
evaluated by using the steepest descent method twice. One finds the following analytical 
estimates 

1 



Y 



Yi2{z>za) ^ lf-^-^KuTt[Y{zA)] 

gi(2-2A)Re(A)g-(2--2A)|Im(A)| ^ 

g-i(2-ZA)Re(A)g-(2-2A)|Im(A)| 



(32) 



(33) 



(34) 



where ~ l/Im(A) and the z^s are implicitly defined by Kaah{za) — 1- We remark that 
the relation Yp^ = (Yajs)* holds, this assures that the diagonal asymmetries are real. To a first 
approximation we can take zi Z2 = z, which is true up to logarithmic corrections. From 
Eqs. (|32p - (j34p it is possible to find an expression for the trace of Y , which allows us to write 
the diagonal asymmetries explicitly: 



Yii 



g-(2-2A)|Im(A)| 



2g,z \ Kii KuK (K11K22 - K12K21 cos [{z - ZA)Re(A)] e-(2-A)|im(A)|) 



{euK22 + e22Ku)Ki2K2i cos [{z - 2A)Re(A)] - — i^iii^22 ( K2ieue'^^-'^'>^<^^ + c.c. 

z\ 



1 f eno g-(2-2A)|Im(A)l 
Yoo — < \ z X 

2g,z \ K22 K22K {K11K22 - K12K21 COS [{z - ZA)Re(A)] e-(2-^A)|im(A)|) 

(eiii^22 + e22i^ii)Js:i2Js:2i COS [{z - ZA)Re(A)] - —K11K22 {K2iei2e'^^-^'^'^^<^^ + c.c. 

ZA ^ 



(35) 



(36) 



The terms proportional to ^aa/K^a are the familiar asymmetries in the strong wash-out regime, 
while the remaining terms are the corrections due to the correlation between flavours. Such 
corrections are quickly damped by the imaginary part of A, and this behaviour is also confirmed 
by numerical simulations. In the limit A — > 00 we recover the total lepton asymmetry of two 



11 



decoupled flavours: 



Tr[y] = yn + Y,, ^- + f-) . (37) 

2g^z \Kii K22J 

as expected. On the other hand, the limit A — > leads to 

Tr[y] — > — ^ -— ; -— ^ . (38) 

2g^z A11A22 - -Ki2ii2i 

It is easy to see that the quantity on the right hand side is left invariant by a transformation 
of the matrices K, e of the form 

K MKN , e ^ MeN , (39) 

where M, are two generic 2x2 non-singular matrices. In fact, the denominator in (j38|) is 
just the determinant of K which simply transforms as: det(i^') det(M) det(A^) det(i^). On 
the other hand, the numerator may be written as eijEmnKim^jn, where e is the antisymmetric 
Levi-Civita symbol in two dimensions and summation over repeated indices is assumed. So, 
the numerator in (1381) transforms as: 



— ^ij^mn{-^iaKab^bm)(,Mjp^pg-^qn) — 
= {£ijMiaMjp){emnNhmNqn)Kab^pq = 

= det(M) det{N)£apebqKabepg 

= det(M) detiN)eijemnKrmejn (40) 

under (j39p . Therefore the numerator picks up an extra factor, namely det(M) det(A^), which 
exactly cancels that in the denominator and the invariance of (j38p is proved. This fact means 
that a transformation of K and e matrices does not affect the trace in (|38|) . In particular, we 
can evaluate it in the rotated flavour basis defined in Section [3l and obtain 

Tr[Y] ~ = ^ , (41) 

KuK22 K 

which is the single-flavour result, as expected. In the one-flavour limit. Mi ^ 10^^ GeV, the 
efficiency factor r]{K) for the final baryon asymmetry depends only upon K. In the opposite 
limit. Ml ^ 10^^ GeV, the final baryon asymmetry depends upon two different efficiency 
factors, one for each K^a- As discussed in [18], Kaa/K < 2 for large mixing angles and 
therefore the efficiency is enhanced by 0(2) when going from Mi 3> 10^^ GeV to Mi <^ 10^^ 
GeV. 
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Figure 1: The ratio between the lepton asymmetries Yn (green) and (red) computed including 
the off-diagonal terms of Eqs. fl28l) and fl29|) and the ones neglecting them (see text) as a function 
of Ml (left). The trace of the lepton asymmetry divided by the same trace computed in the single- 
flavour approximation (see text) as a function of Mi (right). The parameters are K = 50, Kn = 40, 
K22 = 10, K12 = K21 = 20, eii = 0.4, 622 = 0.1, ei2 = £21 = 0.2. Here and in the following, the 
relative magnitudes of the e entries are chosen consistent with the bounds (fTTI) . 

Figure[T]on the left shows Yaa /{yaa) dec, the diagonal lepton asymmetries Yaa, as functions 
of Ml. In this figure, as well as in all others, we have chosen compatible values for the 
parameters K^fs by fixing the Yukawa couplings Ajq. The analytical results reproduce the 
numerical ones within 10%. On the right we show Tr[y]/yi -flavour as a function of Mi. We 
see that the ratio tends to unity for Mi ^ 2 x 10^^ GeV in agreement with our previous 
findings. In our numerical example, the two flavours give rise to the same asymmetries, 
and for the bound discussed in Section S] to be in the full two-flavour state would require 
Ml ^ (i^22/^)10^^ GeV ~ 2 X 10^^ GeV. However, we see from our numerical results that the 
two-flavour state is reached for larger values of Mi. To our understanding this is due to the 
rapidly oscillating behaviour of the off-diagonal terms. As we already mentioned, computing 
the flavour asymmetries involves an integral over time (or, better, over the parameter z). 
Since the quantum correlators not only decay, but also have a rapid oscillatory behaviour, this 
restricts the range of time integration, thus leading to a suppression of the contribution from 
the flavour oscillations. This effect is magnified by the fact that the oscillations have a time 
scales which is about a factor of ten smaller than the damping timescale. We deduce from our 
results that even for values of Mi ~ 10^^ GeV the full two-flavour regime is attained. 
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Figure 3: The time evolution of the asymmetries for Mi 
Ki2 = K21 = 20, en = 0.4, £22 = 0.1, = £21 = 0.2. 
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Figure 4: The time evolution of the asymmetries for Mi = 5 x 10^^ GeV, 
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In Figs, m [3] and m we present the evolution of the asymmetries for a given choice of the 
parameters. As expected, for smaUer values of Mi the off-diagonal terms die out for larger 
values of z. However, by the time the asymmetries stored in the diagonal terms are frozen out, 
the flavour oscillations have already been wiped out. 



5.2 Strong wash-out for one flavour and weak wash-out for the 
other one 

This regime is characterized by K22 ^ 1 ^ -f^ii- The main contribution to the total decay 
parameter comes from the strongly interacting flavour K ~ K\\ S> 1, which means that A^i's 
are almost in equilibrium, as in the previous case. Since the damping of the off-diagonal terms 
is sensitive to it is still possible to perform the integrals for Y\\ and Y22 by means of the 
steepest descent method, getting the same estimates as in the previous regime. We find 
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satisfies the property /(A — > 00) = 0. As in the previous case, one first finds an expression for 
the trace of Y and then uses it to write the diagonal entries in an explicit form 
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Figure 5: The ratio between the lepton asymmetries Yn (green) and (red) computed including 
the off-diagonal terms of Eqs. fl28l) and fl29|) and the ones neglecting them (see text) as a function 
of Ml (left). The trace of the lepton asymmetry divided by the same trace computed in the 
single-flavour approximation (see text) as a function of Mi (right). The parameters are Ku = 30, 
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10-2, Ki2 = K21 = 0.6, en = 0.3, 622 = 5 x lO'^, eu = 621 = 0.006. 



If A ^ 00 the previous expressions reduce to those usually found in the literature [TT] , where 
the off-diagonal correlations are neglected and the two flavours are completely decoupled 



Yii 



1 



en 
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(47) 



Figure [5] on the right shows Tr[y]/yi_fiavour as a function of Mi. We see that the ratio 
tends to unity for large values of Mi, as expected and it does it very fast, in agreement with our 
previous findings that, as soon Ml > 10^2 Qgv, then the two-flavour regime is reached. Figure 
[5] on the left shows Yaa/(Yaa)dec as functions of Mi. The analytical results reproduce the 
numerical ones within 10%. From this figure we deduce that neglecting the off-diagonal terms 
in evaluating the diagonal terms of the matrix y is a good approximation for the strongly 
washed-out flavour for values of Mi ~ 10^^ GeV. For the weakly coupled flavour the transition 
occurs at Mi > {10/K)10^^ GeV ~ 3 x 10^^ GeV, as derived in Section 4. This time the 
transition does not occur for values of Mi ~ 10^^ GeV because, for the set of parameters 
chosen, the asymmetry stored in the weakly coupled flavour is comparable with the one stored 
in the off-diagonal terms. This illustrates the fact that the contribution from the off-diagonal 
terms may influence the final asymmetry in the weakly coupled flavour if the choice of the 
parameters is such that the off-diagonal CP asymmetries and wash out factors are not too 
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Figure 6: The ratio between the lepton asymmetries Yn (green) and (red) computed including 
the off-diagonal terms and the ones neglecting them (see text) as a function of Mi for Ku = 2.4, 
K22 = 0.6, K12 = K21 = 1.2, en = 0.25, £22 = 0.06, = £21 = 0.12. 

small. This might be relevant if the weakly coupled flavour gives the largest contribution to 
the final baryon asymmetry. On the other hand, one would expect that, when the asymmetry 
stored in the weakly coupled flavour is large enough, then the values of ^22 computed with 
and without taking into account the off-diagonal terms should be very close. This expectation 
is shown to be correct in Figure [H It illustrates also our previous estimates that, if K^a ~ 3, 
then the full two flavour regime should be recovered for Mi < (3/K)10^^ GeV ^ 10^^ GeV. 
In Figs, d [8] and [9] we present the evolution of the asymmetries for a given choice of the 
parameters. Again, for large values of Mi the off-diagonal terms die out for larger values of z. 
However, by the time the asymmetries stored in the diagonal terms are frozen out, the flavour 
oscillations have already been wiped out. 
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Figure 8: The time evolution of the asymmetries for Mi = 10 GeV,, K ~ 30, Ku = 30, K22 
10-2, Ku = K21 = 0.6, en = 0.3, e22 = 5 x 10-^, ei2 = eai = 0.006. 
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Figure 9: The time evolution of the asymmetries forMi = 5 x 10^^ GeV, K ~ 30, Ku = 30, 
K22 = 10-2, K12 = K21 = 0.6, en = 0.3, e22 = 5 x 10-^, ei2 = e2i = 0.006. 
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6 Comments and Conclusions 



In this paper we have studied the impact of the oscillations among the lepton asymmetries in 
leptogenesis and investigated the transition from the one- flavour to the two- flavour states. We 
also accounted for the AL = 1 scatterings both in the CP asymmetries and in the wash-out 
terms. The transition mimics the realistic one when the r flavour becomes distinguishable from 
the other two flavours. We have first formally shown that for Mi ^ 10^^ GeV, the quantum 
correlators are relevant to reduce the system of Boltzmann equations to a single equation for 
the total lepton asymmetry. In this regime the one- flavour approximation holds. Subsequently, 
we have shown that in the regime Mi <C 10^^ GeV, the full two-flavour state is recovered thanks 
to the damping of the quantum correlators. We have subsequently solved both analytically and 
numerically the Boltzmann equations for the lepton asymmetries in flavour space. Particular 
attention has been devoted to the case Mi ~ 10^^ GeV where we expected the role played by 
the quantum correlators to be maximal. 

Let us summarize our results. If all flavours are in the weak wash-out regime, the two flavour 
state is reached and the flavour oscillations may be safely neglected if Mi < lO^^ Q^y jf 
flavours are in the strong wash-out regime, we have estimated analytically that the two flavour 
state is reached and the flavour oscillations may be safely neglected if Mi < {Kaa/^)W^^ GeV. 
We point out however that our numerical studies show that the real bound is weaker. The two 
flavour state is reached even for values of Mi close to 10^^ GeV. The flavour oscillations seem 
to efficiently project the lepton state on the flavour basis. To our understanding this is due 
to the short timescale of the flavour oscillations compared to the damping timescale. Flavour 
oscillations decay and have a rapid oscillatory behaviour, thus restricting the range of time 
integration. This suppresses the contribution from the flavour oscillations to all the dynamics, 
rendering the transition easier. 

We conclude that for the strong wash out case it is a good approximation to solve the 
Boltzmann equations just for the asymmetries stored in the lepton doublets. This procedure 
is usually followed in the recent literature regarding the flavoured leptogenesis. Our results 
justify it. 

The same conclusion is obtained if all the flavours are in the so-called mild regime. This 
occurs when the lepton asymmetry is generated only by the low energy CP violating phases 
in the PMNS matrix [13] . 

In the extreme case in which one of the flavour is very weakly coupled and the other is 
strongly coupled, the approximation of neglecting the flavour oscillations is a good one for the 
strongly coupled flavour even for Mi ~ 10^^ GeV. For the weakly coupled flavour neglecting 
the off-diagonal terms may be too drastic for Mi ~ 10^^ GeV, especially if the parameters of 



19 



the off-diagonal terms are such that they induce large asymmetries. However, as soon as Mi is 
smaller than the analytically estimated value ~ (10/-fC)10^^ GeV, neglecting the off-diagonal 
terms is safe. 

Our findings therefore indicate that the flavour effects in leptogenesis become generically 
relevant at Mi 10^^ GeV. Let us conclude with some comments. In this paper we have 
dealt with classical Boltzmann equations. However, a full treatment based on the quantum 
Boltzmann equations would be welcome to study in detail the transition from one- to the 
two-flavour state. A full quantum treatment usually introduces memory effects |19] leading 
to relaxation times which are longer than the one dictated by the thermalization rates of the 
particles in the plasma. In the quantum approach, particle number densities are replaced 
by Green functions. The latter are subject both to exponential decays and to an oscillatory 
behaviour which restrict the range of time integration for the scattering terms, thus leading to 
larger relaxation times and to a decrease of the wash-out rates. This might further help the 
flavour oscillations to efficiently project the lepton state on the flavour basis. 

If the RH spectrum is quasi-degenerate, leptogenesis takes place through a resonance effect. 
In such a case the final baryon asymmetry does not depend any longer on the mass of the RH 
neutrinos. Therefore, Mi may be chosen to well reproduce the full flavour regime without 
causing any suppression in the final baryon asymmetry. 

Finally, let us comment about the upper bound on the neutrino mass from leptogenesis. In 
the one-flavour approximation there is a bound on the largest light neutrino mass m because the 
total CP asymmetry is bounded from above. The upper limit scales like Mi/m [20]. Therefore, 
larger values of m needs larger values of Mi to explain the observed baryon asymmetry. 
However, Mi may not be increased indefinitely, because at Mi ~ (eV/m)^ lO^'' GeV, AL = 2 
scatterings enter in thermal equilibrium and wipe out the asymmetry. This leads to the upper 
bound m ^ 0.15 eV. In flavour leptogenesis the bound on the individual CP asymmetries 
(jlip scales like m and therefore it was concluded that no bound stringent exists on the largest 
light neutrino mass [9j. From these considerations it is clear that the bound on fn depends 
very much on which regime leptogenesis is occuring, i.e either the one-flavor or the two-flavour 
regime. For large values of m, the strong wash-out regime applies and, as we have seen in 
Sec. m the full flavour regime roughly (because our numerical results indicate that the bound 
is weaker) holds only for Mi < (Kaa/ K)10^'^ GeV. Therefore, one would expect that, again, m 
cannot be large at will since K scales as fn. Indeed, at m ~ 2 eV the full flavour regime would 
seem not to apply [17]. To get this estimate it is assumed that both flavours are in the strong 
wash-out regime, have roughly the same CP asymmetries, but that one of the two has a wash- 
out coefficient much smaller than the other, 1 <C Kaa ^ -^/3/3- Under these circumstances the 
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final baryon asymmetry is dominated by the flavour a 



where we have applied the upper bound (jlip and remind the reader about the bound on M\ 
for the full flavour regime to hold. Since the upper bound is inversely proportional to K^a^ the 
most favourable value for the wash-out factor of the flavour a in the strong wash-out regime 
is Kaa ~ 3.3. Therefore, the maximal baryon asymmetry would be 



5* Svru 



2 (49) 



Setting K ~ (m/0.5 x 10^^ eV), we reproduce the statement that for m ^ 2 eV one is entering 
the one-flavour regime ^17]. This conclusion would seem to indicate that a bound on the 
light neutrino mass m from leptogenesis might be present (even though not useful, given the 
conservative upper bound m < 2 eV from cosmology [21]). We notice, however, that upper 
limit on M\ to be in the two-flavour regime becomes weaker if all flavours have the same wash- 
out term. Assume that the total CP asymmetry e\ is very close to zero (for exactly degenerate 
light neutrino masses ei = and eaa = "f/?/?)- As before, all flavours are in the strong wash- 
out regime, but this time we suppose that — [2]. Under these circumstances the final 
baryon asymmetry reads 

where the flavour a can be identified with the r-flavour and we have applied the formulae 
in Ref. [11] which account for the connection among the asymmetries in the lepton doublets 
and the ones in the charges. Taking K^a/K ~ 1/2, Kaa — {m/0.5 x 10^'^ eV), and, for 
instance. Mi ^ 5 x 10^'^ GeV (which is much larger than 10^^(3/K) GeV ~ 10^ GeV), we are 
well in the full flavour regime. Using the condition (jlip . the following maximal value of the 
baryon asymmetry is achieved 



0.16 

_] 

m 



YB^6i^] X 10-" , (51) 



It shows that, even for light neutrino masses in the few eV range, a large baryon asymmetry 
is generated. We therefore conclude that the bound on the largest of light neutrino mass is 
evaded in flavour leptogenesis. 
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